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Finite-size scaling functions are investigated both for the mean-square magnetization fluctuations 
and for the probability distribution of the magnetization in the one- dimensional Ising model. The 
scaling functions are evaluated in the limit of the temperature going to zero (T — > 0), the size of 
the system going to infinity (N — * oo) while N[l — tanh( J/ksT)] is kept finite (J being the nearest 
neighbor coupling). Exact calculations using various boundary conditions (periodic, antiperiodic, 
free, block) demonstrate explicitly how the scaling functions depend on the boundary conditions. 
We also show that the block (small part of a large system) magnetization distribution results are 
identical to those obtained for free boundary conditions. 
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I. INTRODUCTION 

Finite-size scaling has been developed intensively dur- 
ing the last few decades EI E HI 3 an( l it has be- 
come a standard tool in the studies of critical sys- 
tems. An interesting application of the method is using 
the finite-size scaling of the distribution function of the 
order-parameter fluctuations as hallmarks of universal- 
ity classes. The idea goes back to Bruce 0, 0, 13 wn0 
used it, e.g., to verify that the gas-liquid transition of 
the two-dimensional Lennard- Jones fluid belongs to the 
Ising universality class 8]. The list of applications to 
thermodynamic critical points is long H and the idea 
has reemerged in nonequilibrium critical (or effectively 
critical) systems ES El 13 El , 

as well. 

The usefulness of scaling functions as hallmarks of uni- 
versality classes depends on their availability for com- 
parisons. Indeed, a significant portion of the applica- 
tions is about the Ising universality class where the scal- 
ing functions for the d = 2,3,4 dimensional distribution 
functions are well known from simulations |l4l ITU and 
field-theoretic results in d = 4 — e are also available El • 
The picture gallery of scaling functions is clearly far from 
complete since these functions have been systematically 
worked out only for surface growth models |I0| and for 
Gaussian 1/ f a type noise processes El- 

An important issue concerning the critical distribution 
functions is their dependence on the boundary conditions 
(BC). While theoretical calculations usually address pe- 
riodic systems, the building of a histogram of a physical 
quantity (e.g., the magnetization) in an experimental sys- 
tem involves measuring the magnetizations in patches of 
a given size within the bulk of the system (correspond- 
ing to block-spin magnetizations in an Ising model). As 
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has been shown by Binder [Tj], the block-spin distribu- 
tion function at criticality depends on the BC on the 
block, a finding that is not entirely unexpected since the 
infinite-range critical correlations feel the boundaries of 
the system. 

The BC dependence of the scaling functions is an in- 
teresting problem and, indeed, there has been a series of 
works where the PDF of the magnetization for the d = 2 
Ising model at the critical point has been investigated for 
various BC E3 including some exotic ones (Mobius strip, 
Klein bottle). Similar problems have also been studied 
for the roughness distribution of 1/ f a type noise pro- 
cesses El- Analytical results about BC dependence of 
the scaling functions are scarce: they are restricted to 
Gaussian models |l3j| . expansions around d = 4 |l6j| and 
around the spherical limit |l8j |. This is why we decided 
to revisit the the d = I Ising model where the effect of 
BC can be seen in analytical detail. 

Although the critical temperature of the d = 1 Ising 
model is zero, it displays nontrivial features in its finite- 
size scaling as the critical point is approached. The par- 
ticular case of the distribution function of the magnetiza- 
tion in bulk blocks has already been discussed by Bruce 
d- The purpose of this paper is to calculate the scal- 
ing function for the case of periodic (PBC), antiperiodic 
(APBC), free (FBC) and block (BBC) boundary condi- 
tions, and thus gauge the importance of the role played 
by the BC. 

For pedagogical purposes, we also compute the finite- 
size scaling of the magnetization fluctuations. The cal- 
culation is elementary in this case and one can easily ob- 
serve that the periodic, antiperiodic, and free BC yield 
distinct scaling functions. Furthermore, one can also see 
explicitly how the scaling function associated with the 
block BC emerges when a small part of a large system is 
used for measuring the fluctuations. 

The evaluation of the magnetization distribution is 
somewhat more involved but can be carried out relatively 
simply by using a combinatorial approach. For the peri- 
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odic and free BC, the calculations yield nontrivial func- 
tions which are combinations of two delta peaks and a 
continuum background, with the relative weight of the 
delta functions reduced for the case of FBC. The delta 
functions disappear entirely for antiperiodic BC. Finally, 
the combinatorial approach reproduces Bruce's result for 
the BBC and it turns out that the block scaling function 
is identical to that of the FBC case. 



II. MODEL AND NOTATION 

We consider the one-dimensional Ising chain of N spins 
(a - , = ±1, i = 1, ..,N) with ferromagnetic (J > 0) cou- 
pling and with various BC. The interaction energy of a 
given configuration of spins {crj} is given by 



N-l 



(1) 



where J(, c = J, —J, for periodic, antiperiodic and free 
BC respectively. The model is exactly solvable using, 
e.g., the transfer matrix formalism |19|| . and the partition 
functions for periodic BC (upper signs) and for antiperi- 
odic BC (lower signs) are 



z (p),(a) = 2 JV( cosh JV K ± sinh N R) 

while the correlations are given by 

v n ± y N-n 



{a, a. 



l±v 



N 



(2) 



(3) 



where v = tannic with K = J/ksT and, furthermore, 
1 < i, i + n < N. The above quantities are particularly 
simple for free BC 
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(4) 



Note that the correlations only depend on the distance of 
the two spins and not on their particular positions within 
the chain. This is true not only in the periodic case but 
also for APBC and FBC. Also note that the (p), (a), and 
(/) superscripts refer to periodic, antiperiodic and free 
BC respectively throughout the paper. 



III. FINITE-SIZE SCALING OF 
FLUCTUATIONS 

The mean square fluctuations of the total magnetiza- 
tion M — J^iLi a i can be calculated via the spin corre- 
lations as 



(M 2 ) = J2 M 



(5) 



We begin with the discussion of (M 2 ) for the periodic-, 
free-, and antiperiodic BC, leaving the case of block BC 
for a separate subsection. 



A. Periodic-, free-, and antiperiodic BC 

Substituting the expressions for the spin correlations 
(l.'il II) into Eq.lJSJ, one easily finds 
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In the thermodynamic limit (TV 
pressions reduce to 



PBC 

APBC (6) 

FBC . 

oo) all the above ex- 



(M 2 
N 



1 



1-v 



(7) 



and 



One can see that the fluctuations diverge as T — > 0, 
they have a power-law singularity provided t = 1 — v = 
1— tanh J /ksT is used as the control parameter. One can 
also read off the value (7 = 1) of the susceptibility 
exponent [2(1 ]. 

The correlation-length exponent, v, is another expo- 
nent needed in finite-size scaling. It is obtained from the 
spin-spin correlations which decay exponentially in the 
thermodynamic limit. The correlation length, £, defined 
by the exponential decay is independent of the BC 



{Pi&i+n/N^oo 



_ y n _ e -nln(l/t)) 



and diverges for T — > as 



1 



In [1 - (1 - v)} 



(8) 



(9) 



thus providing us with v = 1. As noted by Chen and 
Domb [2lJ , the above definition of the correlation length 
(rather than that found from the second-moment of the 
correlation function) is the appropriate one to discuss 
finite size scaling for this model. 

The scaling form one expects in finite-size scaling in 
dimension d is as follows 



lj 



(M 2 



<$>(tN 



l/dv\ 



(10) 



with N being the number of spins. Putting in the above 
expression d = 7 = v = 1, the finite-size scaling sugges- 
tion takes the form 



(A/ 2 ) 
N 2 



(11) 



where m = M/N is the magnetization density. 

Dividing both sides of Eqs. © by N, we indeed find 
that (M 2 )/N 2 = (m 2 ) yields well defined scaling func- 
tions in the limit of t = 1 — v — > while tN — 2£ is kept 
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finite. The scaling variable £ has a simple meaning since, 
in the scaling limit 



2C = Nt = N/£ 



(12) 



i.e., 2£ is the average number of domains (or domain 
walls) in the system. The actual scaling functions are 
given below 



(M 2 ) 
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PBC 



APBC (13) 



FBC, 



and they are shown on Fig. ^ The functions belonging 
to different BC are clearly distinct thus demonstrating 
explicitly the BC dependence of the scaling functions. 



We close this section with a note on the speed of con- 
vergence of the scaling function. Looking at Eq.©, one 
can see that (M 2 ) /N converges exponentially to its ther- 
modynamic limit for PBC, while the convergence is only 
power law (A^ 1 ) for FBC and APBC. As to the scaling 
function, it can be easily shown that the convergence is 
power law (Af _1 ) even in PBC case. More precisely, the 
correction term is of the form N~ l g{(,) where g(C) < 
with g(( — > 0) = and — > oo) = — 1. Fig. [3 displays 
(to 2 ) calculated from Eq. for the PBC case (note that 
for finite A* and using the scaling variable £, the results 
are meaningful only for £ < N). It should be clear from 
Fig. |2 that the convergence is slow although it may ap- 
pear to be quite good at small £ due to the particular 
form of g(0- 



0.8 




FIG. 1: Magnetization fluctuations in the scaling limit for 
various BC [see Eq. (I13H 1. The scaling variable is ( = Nt/2. 
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FIG. 2: Correction to scaling for magnetization fluctuations 
for periodic BC. First equation in © is used to calculate 
(rn 2 ) = (M 2 )/N 2 for system with N = 2,4, 8, 16, 32, and oo. 
The scaling variable is £ = 7V£/2.The 1/N dependence of the 
correction can be visually observed (note that the distance 
from the N = oo curve is halved as N is doubled) 



The scaling functions coincide for £ — * oo. This is 
understandable since large £ means large number of do- 
main walls which means that the disordered regime is 
approached where the effects of the BC diminish. The 
differences in £ — > limit can be accounted for by the 
number of domain walls in the system. In particular, 
the ground state is completely aligned for PBC and FBC 
systems while it contains an arbitrarily positioned single 
domain wall in case of APBC. As a result 

( m 2)(p).(/) = l , (m 2 )^ = l/3 (14) 

explaining the values of the scaling functions at £ = Q. 

The small £ = behaviour of (m 2 ) can be under- 
stood in terms of the smallest energy excitations above 
the ground states. These excitations are obtained from 
the ground state by adding a pair domain walls in case 
of PBC and APBC, while they consist of a single domain 
wall for FBC. Their effect is a quadratic (linear) decrease 
of (m 2 ) near C = for PBC and APBC (FBC) systems. 



B. Block (window) boundary conditions 

When studying magnetization fluctuations in an ex- 
periment, one usually divides the system into blocks and 
measures the block magnetizations. The corresponding 
theoretical construction in the d — 1 Ising model is to 
consider a block (window) of length I at a position k in 
a system of total length N, and study the fluctuations of 
the block magnetization defined as 

k+l-l 

(m 2 ) = J2 Ww) ■ ( 15 ) 

Provided the full length of the block is within the system 
(i.e., it does not contain the boundary with the coupling 
Jbc), the correlations entering Ea. H15|) and, consequently, 
(M 2 ) does not depend on the location of the block, k, 
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and we can write 



(Ml) =£ + 2^(^-n)(o-iCTn 



(16) 



71=1 



Substituting now the correlations for PBC one read- 
ily obtains an expression which depends both on the 
block size £ and the system size N: 



2\(p) 



1 



2v(l 



.,N-l 



P 



1 + v N i(l - v) 



1 



i 2 (l 



- V) z 

(17) 

where the superscript (p) denotes the periodic BC. Intro- 
ducing the "aspect ratio" b — £/N, and using the scaling 
variable ( = 1(1 — v)/2, the scaling limit N — > oo, I — > oo 
with b and £ finite yields the following scaling function 



2\(p) 



,?\G0 



e 2 



(C,b) 



(18) 



c 



tanh£/fe - 



-2C 1 



-2(l-b)C/b 



2C 



1 + e- 2 ?/ b 



As one can see, the scaling function goes over into the 
FBC case if b — > while it becomes the scaling func- 
tion for the PBC case if b = 1. The function smoothly 
interpolates between the limiting cases in < b < 1. 

If the block is embedded in a FBC system then (Mf)^ 
can be deduced from the observation that (<7i<7i+ n ) is 
independent of the system size (see Eq. 0}. Namely, 
(Mf)W must coincide with (M 2 ) for the FBC with N re- 
placed by £. Thus we have a scaling function l|13fl which 
does not depend on the "aspect ratio" b 



(M t 



(m|)^(C,6) = K)^(0. (19) 



£ 2 



2 \(/)< 



Finally, the fluctuations in a block embedded in an an- 
tiperiodic chain, can also be calculated and the somewhat 
more complicated result has a similar structure as in case 
of embedding in a periodic chain. Namely, changing the 
aspect ratio from b = 1 to b = 0, the result interpo- 
lates between the APBC and the FBC scaling functions 
in Eq.©. 

The common feature of all the above results is that 
the BC become irrelevant in the limit 6^0 where the 
block is much smaller than the system. Furthermore, 
we find that this bulk behaviour coincides with the FBC 
result. Whether this coincidence with the FBC result is a 
general feature of bulk fluctuations is not quite clear and 
should be investigated in more complicated and higher 
dimensional systems. 



IV. MAGNETIZATION DISTRIBUTION 

If at a given temperature, i.e., at a given correlation 
length, the system size goes to infinity then the magne- 
tization distribution goes to a Gaussian around zero due 



to the Central Limit Theorem, and goes eventually to a 
Dirac delta function 



lim lim Pirn) 



5(0) 



(20) 



On the other hand, at any given system size as the tem- 
perature goes to zero, i.e., the correlation length goes to 
infinity, the magnetization goes to either plus or minus 
one, and the distribution becomes 

lim lim P(m) = ~\6(m + l) +S(m- 1)1 . (21) 

TV^ooT^O 2 

The importance of the C = Nt/2 scaling variable is that 
if N — > oo and T — > in such a way that the correlation 
length is always proportional to the system size (i.e., £ 
is a constant l|12|l ) then a nontrivial distribution arises 
even for the d = 1 Ising model. 

The meaning of £ suggests the development of a small 
C (small number g of domain walls) expansion. Thus 
we shall first calculate P g (M), the probability of a given 
magnetization in the presence of g domain walls. Once 
P g (M) is known, the probability of a given magnetization 
P(M) can be obtained by summing up P g (M) 



P{M)=Y J P a {M) 



(22) 



The states with g domain walls are degenerate as their 
energy E g depends only on the number of walls 



E„ 



J(2g - N) 



PBC, APBC 



J(2 5 -7V + 1) FBC 



(23) 



Thus, the calculation P g (M) reduces to counting all pos- 
sible spin configurations Sl g (M) at given g and M values 



Pg(M) 



-EJk B T 



-Qg(M) 



(24) 



Here Z is the partition function corresponding to a given 
BC which, in the scaling limit, becomes 



'2e NK cosh C PBC 
2e JVK sinhC APBC 
2e (JV-i)^ e C FBC. 



(25) 



For g — 0, i.e., if there is no domain wall in the configu- 
rations 



fl (M) = S m ,n + 8m.- n 



(26) 



For g > 0, let Wj be the position of the jth wall (say 
Wj = i when it is just before ai) with j — 1, ..,g, i.e., 
Wj < Wj + i for j = 1, ..,g — 1, and also note that by 
definition W\ > and W g < N. Let us denote the mag- 
netization of the domain between W\ and W2 by M\. 
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Now, it is sufficient to obtain the number of configura- 
tions with the restriction Mi > and denote it as Q+(M) 



n g (M) = n+(M) + nj(~M) 



(27) 



Note that fixed magnetization also means fixed number 
of upspins 7V T = (N + M)/2 and downspins iVj. = (N - 
M)/2. 



For s 7^ 0, using £ = Ne~ 2K leads to the final result for 
the magnetization distribution with a fixed 2s number of 
walls 



P2s{m) 



(C/2) 2s {l~m 2 ) s - 
s\(s — 1)! cosh C 



(34) 



It becomes clear at this point that we are doing a small 
£ expansion and that a fixed 2 s number of domain walls 
belongs to the order 2s of the expansion. 



A. Periodic boundary conditions 

PBC enforce the number of domain walls to be always 
even g = 2s, with s = 0, [N/2], where [ ] stands for the 
integer part. As there is at least one spin after each do- 
main wall, we can imagine those spins as being attached 
to the walls. A typical configuration looks like 



414 (Wit) TT (^ 2 4) 444... TT (w 2s 4) 44 



(28) 



Now counting all the possible configurations is equivalent 
to distributing in all possible ways the — s (not at- 
tached) upspins among the s up domains (those domains 
where the spins are up, i.e., between Wsj—i and W 2 j) 
and independently distributing the N± — s (not attached) 
downspins among the s + 1 down domains (between Wsy 
and Wzj+i, also in front of W\ and behind Wis) 



n+(M) 



iV T - 1 

s- 1 



(29) 



Note that the binomial coefficient (?) = for a < b, 
which reflects the fact that there are no configurations 
with more domain walls than either 27V j or 2JVj_, i.e., 
s < min(iV|, iVjJ. With the formula (|27|l one can easily 
arrive at without restriction on the sign of Mi 



n 2s (M) = 



N, 



1 



Ni-1 
s - 1 



(30) 



In the N — > oo limit, for fix n-\ = N^/N and n± 
N^/N the number of configurations becomes 



n 2s (M) = N 



s!(s-l)! ■ 



(31) 



In this limit we also need to switch from the discrete prob- 
abilities P(M) to the probability density P(m), which 
brings in a factor N/2, and substituting Eq. (j31(l into 
Eq. J2U leads to 



P2s{m) 



-4sK ]\ T 2s 



N 2s (n m 



lS-1 



(32) 



Now the scaling limit can be finally taken using (|25H . 
and for s = one realizes that the expression of Eq. I|26() 
develops singularities 



Po(m) 



1 



2 cosh £ 



[6(m + 1) + 6(m - 1)] 



(33) 




FIG. 3: Magnetization distribution P(m) in the d = 1 Ising 
model with periodic BC in the T — > and N —* oo limit as 
a function of m = M/N with the scaling variable £ = Nt/2 
fixed at values £ = 1, 2, 4, 8, 16 (from bottom at origin). 
The arrows symbolize the singular parts ~ 8(m ± 1) of the 
distributions of Eq. (1371 . 

The magnetization distribution without restrictions on 
the number of walls can be obtained from Eq. I|22|l which 
is valid also in the continuous limit 



P{m) =^P 2s (m) . 



(35) 



s=0 



Using the series form of the modified Bessel functions 



(x/2) 



2k+v 



k=0 



k\{k + v)\ 



(36) 



leads to the final expression for periodic BC 
1 



P( p \m) = 



2 cosh £ 



[5(m + 1) + S(m - 1)] (37) 



C 



2VT — m 2 cosh £ 



1 — iw 



One can easily check that P(m) is normalized, i.e., 
J_ 1 dm P(m) = 1, by changing the integration variable 
7ti to 9 = arcsin m and using the series form of Eq. 13till . 
In the same way the expression for the second moment 
of P(m) in Eq. (jT3|) can also be obtained. 

One can investigate the speed of convergence in N of 
the magnetization distribution P(m) ofEq. I|37|) . Instead 
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of making Monte Carlo simulations on the Ising model 
we calculated numerically the probabilities of possible 
magnetizations for finite chains, based on Eq. Q24|l and 
(|35[) . and multiplied the results by (N + 1)/2 for the sake 
of comparison with the probability density P(m) in the 
scaling limit 



P(m, AT) 



2Z 



s=Q 



(38) 



where Vt 2s {M) is given by Eq. $3%, CP, and K = 
artanh(l — 2(/N). Eq. (|38[) can be easily evaluated nu- 
merically, as the the computation time increases linearly 
with N as opposed to the exponentially long time needed 
to encounter all possible configurations. One observes in 
Fig. H that the convergence is faster for smaller values of 
C, in agreement with Fig. 




FIG. 4: Magnetization distribution P(m) for periodic BC in 
the N — * oo limit [only the regular part of Eq. 13711 is depicted] 
compared to its finite JV = 20 (square), 40 (circle), and 80 
(triangle) forms of Eq. 1381 with the scaling variable being 
£ = 3 (closed symbols), and 6 (open symbols). Observe the 
evolving singularities at m = ±1 for finite systems (closed 
triangles are out of range). 



B. Antiperiodic boundary conditions 

The main difference from the periodic case is that the 
number of domain walls is odd g — 2s — 1, with s = 
1, [(N + l)/2]. As there has to be at least one domain 
wall in each configuration there are no Dirac delta peaks 
at m = ±1 in the distribution. A typical configuration 
looks like 

14 (Wk T) TTT (W 2 1) || ... || (W 2s -i T) TTT . (39) 

The number of configurations with the restriction of 
Mi > is 



A T - 1 
s- 1 



s-1 



(40) 



as we have — s free upspins to distribute into s up do- 
mains and Ni — s + 1 free downspins for s down domains. 
In the continuous limit the number of configurations be- 
comes 



n 2s -i(M) = 2N 



■2.~ - 2 ("T"l) S 1 
(s-1)! 2 ' 



(41) 



and the magnetization distribution with fixed number of 
walls (^IJ reads 



P 2s _i(m) 



(C/ 2)2s-l(l_ m2) 

(s- l)! 2 sinhC 



2\s-l 



(42) 



Summing this expression up over the possible number of 
walls (|22H leads to the final result 



P( ° )(m) = 2lirih7 /0 ( CV/I ~ m2 ) ' (43) 



C. Free boundary conditions 

In case of free BC the number of domain walls can be 
both even g = 2s and odd g = 2s — 1 with s = 1, [N/2]. 
The only difference in counting all possible configura- 
tions, with a given g, M, and the condition Mi > 0, is 
the supplementary restriction that there always has to 
be a down spin before the first wall for obvious reasons, 
which can be visualized as (| W\ f) m t ne example 12811 . 
Now the number of configurations can be easily obtained 



and one observes that in the scaling limit they are 
equal to the corresponding periodic (|31|l or antiperiodic 



gTJ result, i.e., n { 2 J J(M) = nfJ(M) and £l 2s _ 1 (M) = 

f4s-i(-^0- The energies of both the even and odd states 
are greater than the energies of the corresponding peri- 
odic and antiperiodic states (|23|l by J and thus we can 
write the scaling limit of the probability of a given M as 



>/k B T 



-K 



Zip) 



ZU) 



PW(M) 
P (a) (M) 



(45) 



9 = 2s 

g = 2s - 1. 



Using Eas. l|25|l . one finds that the prefactors of the dis- 
tributions depend only on £, thus collecting the contri- 
butions with different number of domain walls 122f) . we 
obtain an expression for P^\m) through P^^ a \m) 

p(/)( m ) = e ~ c [coshC P (p) (m) + sinh C P (a) (m)] , (46) 
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which leads to 

P^(to) 



— [5{m + 1) + 5{m - 1)] (47) 



2V1 - m 2 



h (C\/l-m 2 



1 — m' 



As we shall see in the next section, the above expres- 
sion applies for bulk BC as well. The bulk result was 
obtained previously by Bruce (see equation (3.20)). 
Note also that the relationship established between the 
distributions for different BC Ij46(l is certainly valid for 
the fluctuations of the magnetization l|13[> as well. 

One should note that the coefficient of the singular part 
is larger for the periodic periodic system is 

more likely to be in the completely ordered steady state. 
This is expected due to the lower energy of the states with 
2s — f domain walls than that with 2s domain walls. On 
FigElwe display the nonsingular part of the distributions. 
This nonsingular part also shows the expected sequence 
from PBC resulting in the most ordered state to APBC 
yielding the most disordered state. 
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FIG. 5: Comparison of the regular parts of the magnetization 
distribution P(m) for periodic, antiperiodic, and free BC in 
the scaling limit, with the scaling variable being £ = 1, 4, and 
10 (from bottom at origin). Note that the Dirac deltas at 
m — ±1 of the periodic and the free case are not displayed. 



between a k and a n=k +i 
P{{a k ,...,a n \) = ^jv-i (48) 



where denotes summing over possible values of 

the spins between sites i and j. One can "integrate out" 
the end spin a\ = ±1 in both the numerator and the 
denominator yielding cancelling factors e K + e~ K . This 
can be repeated till the spin a k is reached and then the 
same can be done starting from the other end (<Xjv) of 
the chain. As a result one obtains 



P({<7fc, ...,cr„}) 



{*•}«.,» 



(49) 



This is the free end probability distribution for the spins 
in the segment [k, k + £], thus the magnetization distri- 
bution P(me), with mi = M?/£, is identical to the FBC 
case given by Eq.(gZI), i.e., P^\mp) = P^\m). 

The above derivation does not hold for chains with 
PBC and APBC and one expects that the P(Mi) depends 
on the aspect ratio b — l/N. In the b — > limit, i.e., 
when the window size is relatively small, however, P(AIe) 
becomes independent from the BC imposed on the whole 
chain. More precisely, this happens in the scaling limit 
£,N -> oo, t -> with C = it/2 and b = i/N kept 
constant, a limit where the correlation length £ ~ £ ^ N 
and the effects of the boundaries of the chain can be 
neglected. 



E. Asymptotic regimes 

The small ( expansion is already given by Eq H22J) . 
For sufficiently small £ the distribution P{m) can be ap- 
proximated by the first few terms, e. g., the Dirac delta 
functions (except for the antiperiodic case) plus a con- 
stant probability density for all BC. 

For large values of C a Gaussian approximation can be 
obtained using only the first term of the large x asymp- 
totic of the Bessel functions 



h(x) « Iq(x) 



(50) 



D. Block boundary conditions 

We start with the simplest case, i.e., impose free BC 
on the chain and investigate the magnetization of a finite 
segment (block) of length £. As we shall show the dis- 
tribution of Mi is identical to the free end result for any 
I. In order to see this, consider the Boltzmann weight 
of a configuration of the spins for the segment of spins 



and using also the condition m«l, where P(m) signifi- 
cantly differs from zero, due to the factor exp(£-\/l — to 2 ) 



P(m) 



2 C/2 



(51) 



P(m) eventually evolves to a Dirac delta function 6(m) 
for £ — > oo. 
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V. FINAL REMARKS 

The fluctuations of the magnetization in the one- 
dimensional Ising model have been investigated near its 
critical point in the limit of T — > and N — > oo with the 
average number of domain walls (2Q = Nt) kept constant. 
A simple combinatorial derivation has been presented for 
the magnetization distributions and it was showed that 
in this limit these distributions are nontrivial well defined 
functions of the magnetization and the control parameter 

C- 

The £ dependence of the magnetization distribution 
should be emphasized, as it means that the order param- 
eter distribution of a system can, in general, depend on 
the way the critical point is approached and the thermo- 
dynamic limit is taken. 

We focused our attention on the effect of boundary con- 
ditions, namely we imposed periodic, antiperiodic, and 
free BC on the chain. The magnetization distributions 
are shown to be sensitive to the BC and well distinguish- 
able for all values of (. For antiperiodic BC the distribu- 
tions differ fundamentally from the periodic and free BC 
case (lack of Dirac delta peaks for APBC). 

We also showed that the distribution of the magneti- 
zation of a segment (BBC) of the whole chain with FBC 



coincide with the distribution of the total magnetization, 
independently of the size of the segment. For PBC or 
APBC chains the above statement is true only if the rel- 
ative segment size b goes to zero. 

It is worth mentioning the analogies to a simple ran- 
dom walk. The corresponding quantity is the distribution 
of the width of the walk (mean square deviation of the 
position of the walker), and for that two distinct func- 
tions have been exactly derived for periodic and free BC 
[l3[ . The qualitative behaviour of BBC is also the same 
as that of the Ising model discussed above. 

It would be very interesting to see similar calculations 
for the order parameter distribution in other equilibrium 
one-dimensional models, e.g., the classical and quantum 
XY and Heisenberg models, or nonequilibrium ones, e.g., 
absorbing state phase transitions. 
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